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Abstract 

The current phenomenological determinations of as (Mt) and as {Mz) are shown to be only 
marginally consistent with the QCD evolution of the strong coupling constant between Mz and Mr. 

This motivates a revised estimate of as {Mr) since the perturbative series used to extract as {Mr) 
from the r hadronic width exhibits slow convergence. Fade summation techniques provide an estimate 
of these unknown higher-order effects, leading to the revised determination as {Mr) = 0.333 ± 0.030. 

This value is 10% smaller than current estimates, improving the compatibility of phenomenological 
estimates for as {Mr) and as {Mz) with the QCD evolution of the strong coupling constant. 

The Particle Data Group (PDG) quotes the following values for the strong coupling constant as deter¬ 
mined from and r decays Q. 

as {Mr) = 0.370 ± 0.033 (1) 

a^(Mz) = 0.118 ±0.003 (2) 

Since these determinations of Og occur at such widely separated energies, the compatibility of these 
values of a<j with the QCD evolution of the coupling constant is an important test of both QCD and the 
phenomenological results used to extract the coupling constant from the experimental data. In particular, 
Us {Mr) is sufficiently large that presently unknown terms from higher order perturbation theory could 
substantially alter the value of {Mr) extracted from the experimental data. Pade approximant methods 
provide estimates of the aggregate effect of (presently unknown) terms from higher-order perturbation 
theory p, |5|. As shown below, the use of Pade summation to estimate such terms leads to a 

substantial decrease in the value of Og {Mr) extracted from r decays, improving the compatibility of 
Ug {Mr) and Og {Mz) with the QCD evolution of the coupling constant. 

The QCD evolution of the coupling constant is governed by the (3 function which is now known to 
4-loop order Q. Using the conventions of Q, a = ^ satisfies the differential equation 
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Using the value as (M^) as an initial condition, the coupling constant can be evolved to the desired energy 
using the differential equation (P) . The only subtlety in this approach is the location of ffavour thresholds 
where the number of effective ffavour degrees of freedom Uf change. In general, matching conditions must 
be imposed at these thresholds to relate QCD with Uf quarks to an effective theory with n/ — 1 light 
quarks and a decoupled heavy quark |^]. Using the matching threshold fith defined by mg{nth) = 
where rriq is the running quark mass, the matching condition to three-loop order is 




1 + 0.1528 


^ + {0.9721 - 0.0847 (n/ - 1)} (fith) 


( 6 ) 


leading to a discontinuity of as across the threshold. Thus to determine the coupling constant at energies 

(5) 

between the c quark threshold and the b quark threshold, the /3 function with nj = 5 is used to run as ^ 
from Mz to Hth = rribinth) = rub using (D as an initial condition. The matching condition (|) is then 
imposed to find the value of ai'^^(m;,) which is then used as an initial condition to evolve as to lower 
energies via the rej = 4 /3 function. 

If as {Mz) is used as the input value to determine the QCD prediction of as (Mt), then one might 
legitimately be concerned about the effect of (unknown) higher-order terms in the (3 function at lower 
energies where as is larger. Fade approximations have proven their utility in determining higher-order 
terms in the (3 function. For example, using as input the four-loop (3 function in 0{N) gauge theory 0, 
asymptotic Fade methods described in Section II of Q are able to predict the five-loop term to better 
than 10% of the known five-loop contributions for < 4 i0]- When these same methods are applied 
to QCD, the following predictions for the unknown five-loop contribution to the (3 function are obtained 

0 - 

nf = A: /34 = 83.7563 (7) 

nj = 5: /34 = 134.56 (8) 


From these predictions, f3 functions containing [2|2] Fade approximants can be constructed to estimate the 
sum of all higher-order contributions. These Fad+summations, whose Maclaurin expansions reproduce 
/?!, (32, (33 and and the asymptotic Fad+approximant estimates (0,^) of (3^, are given by: 


n/ = 4: 

/3(a) 
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nf = 5: 

/3(a) 

23x' 

12 


1 - 5.8963a-4.0110a2' 
1 - 7.4363a + 4.3932a2 

1 - 5.9761a - 6.9861a2 
1 - 7.2369a - 0.66390a2 


(9) 

( 10 ) 


Thus the QCD prediction of a^ (Mr) depends on only two parameters: the initial condition a^ (M^) 
and the position of the five-flavour threshold defined by mb{mb) = mb- As will be discussed below, the 
uncertainty in the Farticle Data Group value |Q] for this threshold 

mb{mb) = (4.3 ± 0.3)GeV (11) 

has a negligible effect on the QCD prediction of a* {Mr) compared with the uncertainty in a^ {Mz) (§). 
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The compatibility of the experimentally/phenomenologically determined values a* (Afz) and Ug {Mr) 
with the QCD evolution of the coupling constant can now be studied. Figure shows the effect on as{Q) 
of progressive increases in the number of perturbative terms in the /3 function, culminating with the Fade 
summation (^,^) for (3. It is evident that the curves for as{Q) appear to converge from below to that 
generated by the Fade summation of the {3 function, since the gaps between curves of successive order 
decrease. Using the input values (|,11) for the QCD evolution of down from the to r mass, we 
obtain the following range of values for Ug {Mr) for successive orders of perturbation theory [pi]]: 


2 - loop 0.2910 < {Mr) < 0.3391 

3 - loop 0.2944 < {Mr) < 0.3451 

4 - loop 0.2957 < {Mr) < 0.3477 

Fade summation 0.2963 < ag {Mr) < 0.3489 


( 12 ) 

(13) 

(14) 

(15) 


The dominant effect on the uncertainty in the above originates from ag {Mz )— the effect of the uncer¬ 
tainty in the five-flavour threshold (0) is inconsequential. Thus as illustrated in Figure the empirical 
determinations of ag at Mz and Mr are only marginally consistent with the QCD evolution of the coupling 
constant. This motivates a revised determination of ag {Mr) using Fade approximation techniques. 

The seminal QCD analysis of the r hadronic width lu 
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leads to the following result for Rr. 
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To facilitate comparison of the Fade improvement of this result with previous research we use values for 

(18) 


the parameters in (17) identical to |12|: the CKM matrix elements, 
|Krf|= 0.9753 , |K,|= 0.2210 

and the Cabibbo angle 


sin^ 6r = 


IK, 


\Vudf + \Vusf 


The electroweak contributions in ( pTj) 
Sew = 1-0194 , S'ew = 0-001 


are 


(19) 


( 20 ) 


The mass-independent perturbative contribution to Rr is the dominant effect: in the MS scheme, it 
takes the form 

1 + 5^^^ = l + a{Mr) + 5.2023[a{Mr)f+ 26.366[a{Mr)f . (21) 

Note that for a characteristic value ag {Mr) = 0.37, the contributions of various orders to are 

1 5 ( 0 ) = 1-^0.118-^0.072 -^0.043 (22) 
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illustrating the slow convergence of the perturbation series, and hence the potentially important role of 
presently unknown terms from higher-order perturbation theory. 

( 2 ) 

Perturbative quark mass contributions to Rr occur for D = 2 in the summation. The quantity 
representing the average of the axial and vector contributions to R-r, is 



1 -|- —a {Mr) 


m? {Mr) + m3 {Mr) 


(23) 


where {Mr) denotes the running quark mass (of flavour i) which has implicit dependence on ag {Mr) 
through its running behaviour 
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To facilitate comparison of our results with [^] , we use identical values for the RG invariant quark masses 
rhi- 


rhu = 8.7 MeV 


md = 15.4 MeV 


rhg = 270 MeV 


(26) 


The total effect of these D = 2 mass corrections is small in comparison with since the potentially 
large s quark contributions are suppressed by sin dc in (^) . 

Dimension-four {D = 4) and higher contributions to Rr are dominantly non-perturbative in origin, 
and involve the QCD condensates which originate from the operator-product expansion [^]. The average 
of the axial and vector contributions for D = 4 is 
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(27) 


which has explicit and implicit dependence on ag {Mr). Values for the QCD condensates appearing in 
(^^, as used in [^], are 

{aGG) = 0.02 GeV^ (28) 

(mj'ipj'ipj) =, /Irf =/itj = 189 MeV , /is = 160 MeV (29) 


Dimension-six contributions are the dominant non-perturbative effect. After use of the effective 
phenomenological equivalent of the vacuum saturation hypothesis |T^] parametrized by the (near) RG- 
invariant quantity, 

pagiip^)^ = 3.8 X 10"^GeV® (30) 

the average of the axial and vector D = 6 contributions is 

^ ^ -0.007 (31) 
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Dimension-eight and higher dimensional non-perturbative contributions to Rt are negligible 

Fade approximation techniques can now be applied to the determination of as {Mr) from an experi¬ 
mentally determined Rr- Asymptotic Fade approximation (AFAF) techniques use convergence properties 
of the perturbative expansion to provide improved estimates of the fifth order term in a perturbation 
series given knowledge of the fourth-order series |Q] . Explicit field-theoretical tests of the AFAF approach 
where the fifth order terms are known lead to extremely accurate results. For example, the A^ = 1 limit 


of the five-loop {3 function for 0{N) gauge theory has been calculated explicitly []^|. Using the AFAF 
algorithm described in Section II of Q, it is possible to make an AFAF estimate of the known five-loop 
term (3^ based upon the contributions from the previous four terms f3i-(3i. This estimate has been found 
to be well within 0.2% of the answer obtained by direct calculation |10|. 


Applying the same AFAF algorithm to the perturbation series (^) leads to the following prediction 


for the [a (Mr)] perturbative contribution: 

1 + 5(0) = 1 + a (Mr) + 5.2023 [a {Mr)f + 26.366 [a {Mr)f + 132.44 [a (M^)]' 


(32) 


It is significant to note that this prediction is very close to the maximum estimated size of the fourth 
order effect used to determine the theoretical uncertainty in |T^, indicating an underestimate of the 
higher order effects in previous work. The AFAF result in (^) should be compared with that obtained 
previously from (non-asymptotic) Fade methods [0 


1 5(0) = l + a (Mr) + 5.2023 [a {Mr)f + 26.366 [a {Mr)f + 109.2 [a {Mr)Y 


(33) 


The coefficient of the term is 20% smaller in the non-asymptotic prediction. 

To determine the effect of the AFAF prediction (32) on 5^*^) we again consider a characteristic value 
as (Mr) 


= 0.118 -h 0.072 -h 0.043 -h 0.025 


(34) 


The convergence of this series is still quite slow: the 0{a‘^) contribution is 20% of the leading perturbative 
contribution. Such slow convergence indicates that the further higher-order terms could have a significant 
effect on A Fade summation, in this case a [2|2] Fade approximant, provides a means of estimating 


the total effect of higher order terms in the perturbation series. Using the AFAF determination in (^) 
leads to the following [2|2] Fade approximation for 


^ ^(0) _ 1 - 6.5483a (Mr) 10.5030 [a (Mr)] 

~ 1 - 7.5483a (Mr) -h 12.8514 [a (Mr)] 


(35) 


Effects of the Fade improvement of on the q;<j (Mr) dependence of Rr, and hence the determination 
of as {Mr), can now be investigated in both the truncated form (32) or [2|2] summation form (^5[). Figure 
0 compares the dependence of Rr on as {Mr) using the full D = 2,4, 6 contributions [equations (p 



and electroweak effects (|20| ) for the four different scenarios for 6^^^: the known perturbation series 
the non-asymptotic Fade expression (^), the asymptotic Fade (AFAF) expression (^), and the Fade 
summation (|3^ . It is evident that at the FDG value a {Mr) = 0.370/7r = 0.117, the Fade effects lead to 
a significant increase in Rr- In comparison with the perturbative result, the size of the Rr enhancement 
obtained from the [2|2] Fade summation (|^ ) is roughly twice the enhancement obtained directly from 
, indicating the significance of the higher-order effects estimated in the [2|2] Fade summation. To 


(32 
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facilitate comparison with Table 4 in [jT^ ] and for future determination of as {Mr), Table ^ contains Rr 
for selected values of Ug {Mr) in the four scenarios 

The Fade enhancement of Rr implies that the value of ag {Mr) extracted from the experimental 
measurement of Rr will decrease. This is illustrated in Figure which inverts the relation between Rr 
and ag {Mr). Again, in comparison with the perturbative result, the decrease in ag {Mr) following from 

of 


the Fade summation (|35D is roughly twice the decrease one obtains from the Fade estimate (32 
just the next-order contribution to For the perturbatively-extracted FDG value Ug {Mr) = 0.37, the 
Fade summation would lead to a decrease of approximately 10% to ag {Mr) = 0.33. 

In conclusion. Fade summation methods indicate that the FDG central value ag {Mr) = 0.37 should 
be brought down by approximately 10% to include the estimated effect of higher order QCD corrections. 
Assuming the same relative uncertainty characterizing the FDG estimate, we find that 


ag {Mr) = 0.333 ± 0.030 


(36) 


a result in much closer agreement with the ranges (IT,!^) devolving from ag {Mz) via four-loop and 
Fade-improved four-loop QCD (3 functions. Thus, Fade methods are seen to improve the compatibility 
of ag {Mr) and ag {Mz) with the coupling constant evolution anticipated from QCD. 
Acknowledgements: The authors are grateful for the financial support of the Natural Sciences and 
Engineering Research Council of Canada (NSERC). 


^The agreement of the “perturbative” column of the Table with [O provides a consistency check on our calculations. 


6 







References 


[1] Particle data Group, R.M. Barnett et al Phys. Rev. D54 (1996) 1. 

[2] M.A. Samuel, G. Li, E. Steinfelds, Phys. Rev. D48 (1993) 869; M.A. Samuel, G. Li, E. Steinfelds, 
Phys. Rev. E51 (1995) 3911. 

[3] J. Ellis, M. Karliner, M.A. Samuel, Phys. Lett. B400 (1997) 176. 

[4] J. Ellis, 1. Jack, D.R.T. Jones, M. Karliner, M.A. Samuel, Phys. Rev D57 (1998) 2665. 

[5] M.A. Samuel, J. Ellis, M. Karliner, Phys. Rev. Lett. 74 (1995) 3911. 

[6] T. van Ritbergen, J.A.M. Vermaseren, S.A. Larin, Phys. Lett. B400 (1997) 379. 

[7] K.G. Chetyrkin, Phys. Lett. B404 (1997) 161. 

[8] G. Rodrigo, A. Santamaria, Phys. Lett. B313 (1993) 441. 

[9] K.G. Chetyrkin, B.A. Kniehl, M. Steinhauser, |hep-ph/9708255 v2 

[10] H. Kleinart, J. Neu. V. Schulte-Prohlinde, K.G. Chetyrkin, S.A. Larin, Phys. Lett. B272 (1991) 36, 
Phys. Lett. B319 (1993) 545 Erratum. 


[11] V. Elias, T.G. Steele, E. Chistie, R. Migneron, K. Sprague, |hep-ph/9806324 . 

[12] E. Braaten, S. Narison, A. Rich, Nucl. Phys. B373 (1992) 581. 

[13] M.A. Shifman, A.I. Vainshtein, V.I. Zakharov, Nucl. Phys. B147 (1979) 385. 


7 






(Mr) 

Perturbative 

Rr 

Nonasymptotic Fade 

APAP 

Fade Summation 

.20 

3.295 

3.300 

3.301 

3.305 

.21 

3.314 

3.320 

3.322 

3.326 

.22 

3.333 

3.341 

3.343 

3.348 

.23 

3.353 

3.362 

3.364 

3.371 

.24 

3.374 

3.385 

3.387 

3.396 

.25 

3.395 

3.408 

3.411 

3.422 

.26 

3.417 

3.433 

3.436 

3.449 

.27 

3.440 

3.458 

3.462 

3.478 

.28 

3.463 

3.484 

3.489 

3.509 

.29 

3.487 

3.511 

3.517 

3.542 

.30 

3.512 

3.540 

3.546 

3.576 

.31 

3.538 

3.570 

3.576 

3.613 

.32 

3.564 

3.600 

3.608 

3.653 

.33 

3.592 

3.632 

3.641 

3.695 

.34 

3.620 

3.666 

3.675 

3.740 

.35 

3.649 

3.700 

3.711 

3.789 

.36 

3.678 

3.736 

3.748 

3.841 

.37 

3.709 

3.773 

3.787 

3.897 

.38 

3.741 

3.812 

3.827 

3.958 

.39 

3.773 

3.852 

3.869 

4.024 

.40 

3.806 

3.894 

3.913 

4.097 

.41 

3.841 

3.937 

3.958 

4.175 

.42 

3.876 

3.982 

4.005 

4.262 

.43 

3.912 

4.029 

4.054 

4.357 

.44 

3.949 

4.077 

4.105 

4.469 


Table 1: Values of Rt for selected as (Mr). The columns differ only in the treatment of used in 
the “perturbative” column uses (^), the “Fade summation” column uses (^), the “APAP” column 
and the “Nonasymptotic Fade” uses 







































Figure 1: Effect of increasing the order of perturbation theory in the QCD evolution of the strong 
coupling constant using a* (Mz) as an initial condition. Higher-loop terms in the /3 function progressively 
increase a* from the 2-loop order bottom (solid) curve to the Fade summation top (dashed-dotted) curve, 
sandwiching the three- and four-loop curves. 
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Figure 2: Comparison of a* {Mr) extracted from experiment (right data point) and ag (Mr) determined 
by QCD evolution of ag from Mz to Mj- (left data point). 
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Figure 3; as a function of Ug {Mr) using the four different treatments 
the perturbative expression (ID, the dotted curve uses the non-asymptotic 
curve uses the APAP expression (32), and the dashed curve uses the Pade 


of(5(°). The solid curve uses 
Pade (1^ , the dashed-dotted 
summation (^). 
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Figure 4: as (Mr) as a function 
the perturbative expression (ID, 
curve uses the APAP expression 


of Rr using the four different treatments 
the dotted curve uses the non-asymptotic 
(32), and the dashed curve uses the Pade 


of(5(°). The solid curve uses 
Pade (1^ , the dashed-dotted 
summation (^). 
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